Magnetic field effects in energy relaxation mediated by Kondo impurities 
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We study the energy distribution function of quasiparticles in voltage biased mesoscopic wires in 
presence of magnetic impurities and applied magnetic field. The system is described by a Boltzmann 
equation where the collision integral is determined by coupling to spin 1/2 impurities. We derive 
an effective coupling to a dissipative spin system which is valid well above Kondo temperature in 
equilibrium or for sufficiently smeared distribution functions in non-equilibrium. For low magnetic 
field an enhancement of energy relaxation is found whereas for larger magnetic fields the e nergy 
relaxation dec reases again meeting qualitatively the experimental findings by Anthore et al. ( :ond- 



mat/0109297 ). This gives a strong indication that magnetic impurities are in fact responsible for the 
enhanced energy relaxation in copper wires. The quantitative comparison, however, shows strong 
deviations for energy relaxation with small energy transfer whereas the large energy transfer regime 
is in agreement with our findings. 
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I. MOTIVATION AND OVERVIEW 

Energy relaxation of "hot" electrons in disordered con- 
ductors at low enough energies was for a long time be- 
lieved to be determined by direct interaction between 
electrons. Recent experiments on mesoscopic wires (see 
Rcf. |l| and references therein) have shown that Kondo 
impurities (localized spins) lead to much higher energy 
relaxation rates than predicted by the standard theory.H 

On the other hand the inelastic collisions of electrons 
and their spin-flips are directly related to the phase co- 
herence time probed by weak localization effects, such as 
low field magnetoresistance. The problem of decoherence 
in weak loC|ajization was recently revisited and intensively 
discussed.™™ „ _ 

Theoretical studies by various groupstTEj lead to a 
satisfactory and consistent explanation of energy relax- 
ation experiments by Kondo impurity mediated electro 
electron interaction in the gold wires by Pierre et al 
Those gold wires were contaminated by iron impurities, 
and the concentration of the impurities could be inde- 
pendently estimated from the magnetoresistance as well 
as from the temperature dependence of the resistivity. 
As a result, it was possible to carry Qut,a parameter free 
comparison of theory an experiment Dij 

At the same time the copper samples in-the first en- 
ergy relaxation experiment by Pothier et aZJiil were fitted 
using the concentration and Kondo temperature of the 
paramagnetic impurities as free parameters. The param- 
eters obtained from energy relaxation disagree with those 
obtained from-thc magnetoresistance experiments of the 
same sample. E3tj Therefore, if magnetic impurities are 
responsible for these effects is dubious. 
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Since the behavior of magnetic impurities is sensitive 
to the applied magnetic field, studies of energy relax- 
ation in presence of the magnetic field could either rule 
out or validate magnetic impurities as relevant scattering 
process in energy relaxation. Recently, Anthore et al.tj 
reported results of such experiments in Cu wires that in- 
dicate a strong dependence of the energy relaxation on 
the magnetic field suggesting that magnetic impurities 
indeed play a role for the copper wires as well. 

In this article we perform a theoretical study of trans- 
port and energy relaxation in a mesoscopic wire in depen- 
dence on an applied magnetic field. We use a diffusive 
Boltzmann equation to account for the static scatterers 
and focus in the inelastic collision integral on magnetic 
impurities. The findings are in qualitative agreement 
with the experimental data in Ref. [l4] supporting the pre- 
sumption that scattering with magnetic impurities is the 
essential mechanism of energy transfer at low tempera- 
tures. However, the apparent inconsistency between the 
values of the experimentally observed energy relaxation 
rate and the dephasing rate extracted from the magne- 
toresistance in Cu wires remains puzzling. 

Starting with a brief discussion of the experiment we 
propose in Sec. || a simple physical picture to explain the 
anomalous dependence of the energy relaxation on the 
magnetic field. In Sec. Ill follows a theoretical descrip- 
tion in terms of a renormalized Hamiltonian restricting 
the interaction processes to the coupling to electron-hole 



IV the numerical 



pairs only. We then present in Sec. 
procedure and the comparison with experimental results. 
Sec. |v| is devoted to a discussion of the interpretation and 
validity of the approach as well as its possible extensions. 
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II. EXPERIMENTAL SITUATION AND 
PHYSICAL PICTURE 

Here, we briefly describe the experimental situation in 
Ref. |lj and a possible qualitative explanation of their 
findings. 

The experimental setup in Ref. |lj consists of a thin 
copper wire of about 45 nm thickness, 105 nm width, 
and 5 /im length connected with two metallic leads. The 
leads are biased by an external voltage source U — 0.1 
mV and U — 0.3 mV imposing a steady state current 
through the wire. The setup is placed in a dilution re- 
frigerator with a temperature of 25 mK and a magnetic 
field up to 2.2 T is applied. The elastic mean free path 
can be estimated to be much smaller than the length of 
the wire, L, so that the transport of electrons between 
the contacts is diffusive. The diffusion constant D = 90 
cm 2 /s, estimated from the low temperature resistance, 
leads to a diffusion time of td = L 2 /D = 2.8 ns. 

The aim of the experiment was to study energy dis- 
tribution of electrons. The distribution function was de- 
termined by tunneling to an underlying aluminum probe 
electrode. The differential tunneling conductance is given 
by a convolution over the electron distribution functions 
in the wire and in the probe electrode, both for B = 
and B ^ 0, cf. with Refs. |ll] , |l5| and [u], respectively. In 
the case of vanishing magnetic field the aluminum probe 
electrode is in the superconducting state, and the peaked 
density of states of the probe electrode allows one to 
straightforwardly extract the energy distribution in the 
wire from the I-V characteristic of the junction. 

For finite magnetic fields, however, the aluminum 
probe electrode is in the normal state. In this case the 
deconvolution has been made using a zero bias anomaly. 
Unfortunately, the latter procedure is less accurate. Con- 
sidering the numerical transformation depicted in Fig. 2 
of Ref. [l4|, we expect that the strongest variations due to 
uncertainties in the shape and the depth of the zero bias 
anomaly arise at the "Fermi points" e = ±eU/2. Fur- 
thermore, the tunnel probe experiment in equilibrium, 
U — > 0, gave the temperature of 65 mK which differs from 
the actual refrigerator temperature. This might be due 
to an oversimplification of the environmental impedance 
responsible for the zero bias anomaly.L 4 ] Consequently, 
the resulting distribution function has to be taken with 
some care, in particular near the Fermi points. 

The electron distribution function in the absence of 
inelastic scattering is just a linear combination of the 
distributions in the left and right electrodeEj 

/o(e, x) = (1 - x)f F (e - eU/2) + xf F (e + eU/2) . (1) 

Here, /f(c) is the Fermi function while x is the longitu- 
dinal coordinate of the observation point in units of the 
total length, L. At low temperatures the distribution (|l|) 
has steps in both e and x dependencies. These steps are 
smeared by inelastic processes, such as electron-phonon 



and electron-electron interaction. As a result, the distri- 
bution in the middle of the wire turns, put to be almost 
insensitive to the bath temperature. tlt3 The smearing 
depends on the effective inelastic relaxation time, and 
the latter can be estimated from the experimentally ob- 
served distribution function. 

The first experimental 1 ! tB for B = have clearly shown 
that the smearing is too strong to be attributed to the 
electron-electron or electron-phonon interaction. There- 
fore, in the following we do not take these interactions 
into account. ,— ■ 

The qualitative outcome of the experimental is that 
the behavior of the inelastic relaxation rate is a non- 
monotonous function of the magnetic field B. For B < 
B\ = eU/4:.3fj,B, the relaxation rate increases with mag- 
netic field and reaches a maximum at B » B\. At 
stronger fields, B > B\, it decreases with further in- 
creasing magnetic field, and at B « B2 — eU/2/j.B it 
reaches almost the same value as it had at B = and 
then decreases further. 

The explanation of such a complex behavior given in 
Ref. [l4| is based on electron scattering by magnetic im- 
purities. For vanishing magnetic field the spin system 
is degenerate and only second or higher order scatter- 
ing processes contribute to energy relaxation. For finite 
magnetic fields, there exists also a first order contribution 
with energy Eh = 9^bB (equal to the Zeeman splitting) 
transferred to or from the spin system. Therefore, the 
energy relaxation rate increases. However, for Eh > eU 
the spins are completely polarized and can no longer con- 
tribute to energy relaxation. Consequently, again only 
higher order processes are effective. Comparing this ex- 
planation with the experiment, one estimates the gyro- 
magnetic factor for the impurity spins as 3 ~ 2. 

Thus, electron-spin interaction taken into account in 
the lowest order of the perturbation theory explains, 
in principle, the main experimental features. However, 
from the theoretical point of view, there appears a sub- 
tlety. The problem is that the higher order terms, 
estimated within the framework of the i-matrix ap- 
proach in Refs. ||,|7],|l0|, lead to a divergent contribution 
~ J 4 /(e^fEn) 2 to be integrated over. Here J is the renor- 
malized coupling constant which defines the strensth. of 
the electron-impurity interaction. The siiggestionaa'H to 
introduce a cutoff at the Korringa widtho, K ~ J 2 gives 
a result that is comparable with the first order contribu- 
tion. Consequently, there is no systematic expansion in 
powers of the coupling constant, and one needs a gener- 
alized approach which is not based on expansion in terms 
of the interaction strength. 

The aim of the present paper is to develop an ap- 
proach able to treat lowest and higher order contribu- 
tions within a unified scheme. We show that the prob- 
lem can be described as electrons coupled to a dissipa- 
tive spin system. The dissipation of impurity spins is, 
in turn, caused by creation/annihilation of electron-hole 
pairs due to electron-spin coupling. We show that this 
scheme includes the divergent higher-order contributions 
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appearing in the t -matrix approach. Using our approach, 
we derive an electron-spin collision integral expressed 
through spin-spin correlation functions. An important 
feature of these correlation functions is that their depen- 
dences on the electron energy is automatically broadened 
by Korringa-type processes. These processes, however, 
depend on the actual electron distribution rather than 
on the thermal equilibrium distribution as in Ref. [19[ As 
a result, the applied voltage plays the role of an effective 
temperature. 

The crucial difference between the cases B = and 
B ^ is the following. At zero magnetic field all spin- 
spin correlation functions are centered at zero energy and 
behave as K/(e 2 + K 2 ). It is important that in the col- 
lision integral these correlators are multiplied by a com- 
bination of the electron distribution functions which at 
small energies is proportional to e. Consequently, the 
broadening turns out to be unimportantQ, and one can 
omit K 2 in the denominator. This way the resulting colli- 
sion integraLhecomes proportional to J 4 and one recovers 
the resultsQ'Ll of the i-matrix approach. 

At finite magnetic fields, however, the spin-spin corre- 
lation function decomposes into three contributions. The 
non spin-flip part is still centered at zero energy, while 
the two spin-flip contributions are peaked at e = ±Eh- 
Let us for a moment accept the above simplified form of 
the correlation function and assume K — > 0. Since in this 
case K/[{e^E H ) 2 + K 2 } — > Tr8{t^E H ), the two spin-flip 
correlation functions with finite energy transfer lead in- 
deed to the desired first order in J 2 contribution. In this 
simplified case only the non spin-flip correlation function 
contributes to the order J 4 in the collision integral. 

Using the t-matrix approach with a finite cutoff K, in- 
troduced by hand, would lead to a double counting of the 
first order in J 2 contribution. Here, on the other hand, 
the cutoff is included automatically and the first order is 
accounted for correctly. 

In this simplified case the reasoning to explain the ex- 
perimental data follows the lines of the experimental- 
ists. The only difference is that for Eh > eU where 
the spin-flip contribution is already frozen out there re- 
mains just 1/3 of the B = energy relaxation because 
only the non spin-flip component contributes to energy 
relaxation in order J 4 . In practice, K is not constant but 
depends on frequency that the correlation functions are 
not Lorentzian shaped and therefore the spin-flip terms 
will also contribute to the J 4 term. Further, the width 
decreases with increasing magnetic field for the non spin- 
flip component and therefore energy relaxation mediated 
by magnetic impurities dies out for Eh 3> eU. 

We believe that our approach provides a consistent ex- 
planation of the magnetic field dependence of the non- 
equilibrium electron distribution in diffusive wires with 
magnetic impurities. 



III. THEORETICAL DESCRIPTION 

Here we present a simplified version of the theory with 
isotropically renormalized coupling constant J indepen- 
dent of energy. In general, the renormalization can be 
anisotropic and energy dependent. These generalizations 
which do not alter the underlying physics are discussed 
in the Appendix. 

We assume that the metallic wire is in the diffusive 
limit, i.e. the elastic relaxation time is much smaller than 
other time scales. We also assume that the distribution 
function of electrons does not depend on the spin. The 
energy distribution of the electrons is governed by the 
Boltzmann equation 



df(e,x) 1 d 2 f(e,x) 



dt td dx 2 
/{/}= fdu{f(e)[l-f(e-u)]W(u) 
-{l-f(e)]f(e-Lo)W(-Lo)}. 



(2) 



(3) 



Here, we include the density of states p in the scattering 
rate W{lo) and omit for convenience the explicit spatial 
dependence of the distribution function. The rate W de- 
scribes the transitions between two electron states with 
energies e and e — u> mediated by coupling to the dissipa- 
tive spin system. Its explicit form is given by 



W{u>) = (cw p /M) (pJ/2) 2 C(uj) 



(4) 



where c; mp is the impurity density. Further, C(lu) is 
the Fourier transform of a spin-spin correlation function. 
The latter can be split as 



C(t) = [C+(t) + C-(t)]/2 + C z (t) 



(5) 



where 



C ± (t) = {S ± (t)S*(0)), C z (t) = (S z (t)S z (0)). (6) 

The averages here mean the spin and electron trace 
weighted with the unknown non-equilibrium density. 
The time evolution is governed by the Hamiltonian H = 
Hq + Hi where 



kcr C ka C k<r ~ E hS Z 



(7) 



ka 



describes free electrons. Here, operators Ct and C ka 
create and annihilate an electron in a given orbital, fc, 
and spin, cr, state. Cka- is the energy of this state. The 
second term in Eq. (Q) describes a spin 1/2 impurity with 
Zeeman splitting Eh — g^sB. The interaction Hamilto- 



H I = J S-s^C^fffe, 



(8) 



couples electrons to the impurity spin system via the 
renormalized coupling strength, J, rather than the bare 
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one, Jo- Further, the electron is coupled only to one im- 
purity spin since we assume that the impurity density 
Cimp is small enough to neglect higher order terms. 

Using this renormalized Hamiltonian we will restrict 
our calculation for the time evolution of C(t) to cou- 
pling to simple electron-hole pair excitations only, see 
Appendix |X] for details. A similar procedure was already 
used in Ref. |20| to discuss the impurity spin resonance 
linewidth in equilibrium using Baym and Kadanoff 's ki- 
netic equations. Here, the spin-spin correlation functions 
are calculated using the conventional projection operator 
techniqucO. The results read 



C z (u) = 

OtM = 



1 

2u 2 + v z {lu) 2 

2V ±V± {uj) 
[to =F E H } 2 + v±{u 



(9) 
(10) 



The functions v Z) v± describing damping of spin fluctua- 
tions can be expressed through an auxiliary function 



CH= / de'f(e')[l-f(u + e')} 



(11) 



stemming from the coupling to electron-hole pairs. They 
read: 

v z {uo) = 7r(pJ) 2 [P+C(w - E H ) + V_({uj + E H )\ (12) 



v±(u}) = (n/4)(pj) 2 [2C(w T E H ) + CM/Pj 



(13) 



Further, V± is the occupation probability for impurity 
spin up or down, respectively. These probabilities are 
determined by a master equation 

dV±/dt = -T±P± + T T V T , V+ + V- = 1 (14) 

which can be solved in the steady state leading to 

p±-r T /(r + + r_). (15) 

Here, the inverse lifetime for the spin up/down state, T±, 
is determined by the expression 



r, = / ''Ai-'if'iU'] 



(16) 



For a given electron distribution, the set of Eqs. jl5| ) 
and (|l6|) determines the occupation probabilities of the 
spin system. From Eqs. (||) and (||) one can prove the 
sum rule for the correlation function, 

C(t = 0) = J {duj/2n) C{uj) = 3/4 = S{S + 1) , (17) 

which is independent of magnetic field. In the weak cou- 
pling limit 

C±(u) = 2ttT±5{uj tE h ) , C z {lo) = ttS{lo)/2 . (18) 

Inserting this results into the rates (||) and (|T^), we 
recover the Fermi's Golden Rule expressions for the elec- 
tron collision operator for the case of interaction with 



a single spin. The (^-functions here signalize energy con- 
servation. The dissipation leads to an energy uncertainty 
that, in turn, results in a broadening of the ^-functions. 
The final |£csult differs from that obtained by the i-matrix 
approachca by allowing for a finite energy uncertainty. 




0.0 
E/eU 

FIG. 1. Comparison of expe 
bution functions from Ref. [L4| (symbols) with the theoretical 
predictions for various magnetic fields and U — 0.1 mV. The 
thick lines are data gained from solving the Boltzmann equa- 
tion where the impurity spins are out of equilibrium whereas 
the thin ones (differing from the thick lines for two values of 
B only) are determined with impurity spins fixed to equilib- 
rium. The distribution functions are given from bottom to 
top for B = 0.0,0.2,0.4,0.8,1.2 T and shifted vertically by 
steps of 0.2 and 0.3, respectively. 



For vanishing magnetic field, Eh = 0, the occupa- 
tion numbers V± = 1/2, and the correlation function 
simply reads C(u>) — 3C z (oj). Inserted in the collision 
integral, our expressio n w ith the renormalized coupling 
constant J from Eq. ( jB6| ) is consistent with the results 
of Refs. |6|J^. The advantage here is that the cutoff, 
^z(0) = 7r(pJ) 2 £(0), suggested in Refs. is naturally 
included. This cutoff equals to the Korringa width K. 
For weak coupling only elastic scattering survives and we 
get W(e) — r sf ^ 1 (5(e)r-where the time r s f is usually referred 
to as spin-flip timeE3. In this case the collision integral 
vanishes identically and one has to go beyond the lowest 
order. 

Further, our expression for the correlation function in 
equilibrium and vanishing magnetic field coincides with 
the results of Woger and ZittartzcJ based on a Nagaoka- 
like decoupling scheme. Using Suhl's i-matrix for the 
renormalized coupling constants, the spin susceptibilities 
follow from proper kinetic equations. The procedure is 
actually similar to the one used in Ref. ^0|for the analysis 
of the impurity spin resonance linewidth. 
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IV. NUMERICAL PROCEDURE AND 
COMPARISON WITH EXPERIMENT 

For the numerical procedure we use anisotropic and 
energy dependent coupling constants J±(e) and J ± (e) 
given in the Appendix Eqs. (£>£) and (B7), respec- 
tively. They have to be calculated self-consistently as 
functionals of the final non-equilibrium distribution func- 
tion. This procedure leads to a slight complication of the 
formulas in the previous section but does not alter their 
structure. The detailed changes to be made are listed in 
the Appendix ^. Further, in Eqs. (|) and © we fix 
the lifetimes in the denominator at resonance, f z (0) and 
u±(±Eh), respectively, and neglect w-dependence of v z 
and v±. This approximation formally violates the sum 
rule However, this violation appears in higher or- 

ders in J and thus is not worrisome. On the other hand 
in our approximation the Boltzmann equation (||) leads 
to the effective Fermi distribution in the middle of a suf- 
ficiently long wire, i.e., gives the correct "hot electron" 
limit. 



U=0.3mV 




0.0 

E/eU 

FIG. 2. Comparison of experimentally determined distri- 
bution functions from Ref. [14l (symbols) with the theoretical 
predictions for various magnetic fields and U = 0.3 mV. The 
thick lines are data gained from solving the Boltzmann equa- 
tion where the impurity spins are out of equilibrium whereas 
the thin ones are determined with impurity spins fixed to equi- 
librium. The distribution functions are given from bottom to 
top for B = 0.0,0.4,0.8,1.2,1.6,2.2 T and shifted vertically 
by steps of 0.2. 



We have started with the solution ([!]) of the impurity- 
free problem inserted i n Eq s. (B6) and (|B^ ) for the cou- 
pling constants, Eqs. (D5) and (|D6|)_ for the auxiliary 
functions Cz( w ) an d C±( UJ ): an d Eq. (g) for the collision 



integral. Then the distribution function was evolved iter- 
atively using the Boltzmann equation (|2|) with collision 



operator (||). At each iteration both the coupling con- 
stants and the correlation functions were updated. Af- 
ter about 120 iterations a stationary solution has been 
reached. 

To make an independent comparison of the finite field 
data, we have fitted the impurity density with data at 
B = 0. The resulting impurity concentration Ci mp = 
8 ppm is in accordance with the experimental purity of 
copper of 99.999%.E-3 The density of states is chosen to 
be p = 0.21/(site-eV).Ej The Kondo temperature, not 
known so far, has been assumed as Tk = 0.4 K. Further, 
the gyromagnetic factor was chosen as g w 2, see review 
of the experimental results in Sec. ||. 

The comparison of the distribution functions for U = 
0.1 mV is shown in Fig. |l| and for U — 0.3mV in Fig. ||. 
The symbols are the distribution functions for several 
magnetic fields obtained from experimental data by a 
deconvolution procedure. Note, that a magnetic field of 
B = 1 T leads to a Zeeman splitting Eh ~ 0.12 meV. 
Thick lines are the outcome of our numerical procedure 
allowing for the non-equilibrium occupation numbers 
( p"5| ) for the impurity spins. For comparison, thin broken 
lines are the results obtained from the Boltzmann equa- 
tion for equilibrium impurity spins. It is clear that the re- 
sults for relatively weak magnetic fields, Eh < ell, agree 
with theoretical predictions only if the non-equilibrium 
spin population is taken into account. Consequently, we 
conclude that the impurity states are indeed out of equi- 
librium. Since the only spin relaxation mechanism taken 
onto account is an interaction with electrons it follows 
that "hot" impurity centers serve as mediators for the 
electron-electron interaction. 

In large magnetic fields, Eh > ell, according to the 
theory, the impurity-induced energy relaxation is frozen 
out. Consequently, to obtain a quantitative agreement 
with the experimental results shown in Fig. |l| one should 
take into account other, though weak, scattering mecha- 
nisms. 

Further, we observe that the numerical data for U = 
0.3mV and B = 1.2T show an additional step in the 
middle of the energy distribution function. This can be 
explained in terms of lowest order, J 2 , scattering with im- 
purity spins where the transferred energy between elec- 
tron and impurity spin is always ±Eh- The origin is 
that the Zeeman splitting, Eh, at this magnetic field is 
almost ell/ 2. Electrons in the energy region slightly be- 
low —ell/ 2 are scattered slightly below e = because the 
distribution function differs strongly in these two regions. 
The same statement is correct for electrons which scatter 
from slightly above e — to slightly above elf/2. Scat- 
tering between other regions is suppressed by the small 
difference in the distribution function at the contribut- 
ing regions. Using only the first order in J 2 description 
for energy relaxation this feature would be much more 
pronounced as the one shown in Fig. |^. The finite width 
in the correlation functions ( |l0| ) responsible for these in- 
elastic processes smear already out most of the sharp 
features. We assume that additional energy relaxation 
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processes, such as direct electron-electron interaction are 
responsible for the missing step in the corresponding ex- 
perimental data. This feature is therefore assumed to 
be an artifact of our numerical calculation caused by ne- 
glecting other scattering mechanisms. 



U=0.3mV 




FIG. 3. Averaged negative slope of the energy distribution 
functions at the plateau near e as in dependence of the mag- 
netic field B in units of T. The data for U = O.lmV (dots) 
are extracted from Fig. [j] and for U = 0.3mV (squares) from 
Fig. The solid lines are our numerical data and the dashed 
lines the experimental ones. 



Since in our theoretical approach additional scatter- 
ing mechanisms which lead to small energy transfer are 
missing and further the experimental data uncertainties 
are most pronounced at the "Fermi points" e = ±eU/2 
it is not reasonable to consider these energy regions in 
more detail. To characterize the quality of our fits we 
focus the comparison between our numerical data and 
the experiment to the "plateau" region around e w in 
the middle of the two-step distribution function. The 
negative slope of the distribution function is a good in- 
dicator of the energy relaxation strength as long as it is 
small. For vanishing energy relaxation it is zero and it 
increases with increasing energy relaxation. In Fig. || we 
show the averaged negative slope of the energy distribu- 
tion functions at the plateau near e rs in dependence of 
the magnetic field B. We find that the theoretical data 
(solid lines) meet the qualitative outcome in Sec. |l| of 
increasing energy relaxation up to the maximum value 
at about Eh ~ eU/2 and then decreasing again show- 
ing at Eh ~ eU about the same value as at Eh = 0. 
The data, however, lack to explain the experimental out- 
come (dashed lines) in quantitative detail, which is not 
further surprising taking into account the simplifications 
we made. 



V. DISCUSSION 

Let us first discuss the outcome of our findings, then 
focus on the justification of various assumptions not ad- 
dressed in previous sections, and conclude with an out- 
look to further work. 

We have succeeded to explain the main features of the 
inelastic relaxation due to impurity-mediated electron- 
electron interaction - effective relaxation in weak fields 
and its decrease in large magnetic fields, Eh > eU. Fur- 
ther, we found that the impurity spins are out of equilib- 
rium. Namely, for Eh < eU both spin states are occu- 
pied even at zero temperature, and the role of the effec- 
tive temperature is played by the applied voltage. In this 
sense the impurity centers turn out to be "hot" . Contrar- 
ily, when the magnetic field is strong, Eh > eU, only the 
lowest spin level should be occupied, the occupation of 
the "unfavorite" spin state is determined by high-order 
processes. Quantitative agreement in this region, how- 
ever, lacks for small energy transfer, responsible for the 
smearing at the "Fermi points" e = ±eU/2, where the 
experimentally determined distribution functions show 
much stronger energy relaxation than the one expected 
from the theory. It is so far not clear if the small energy 
transfer discrepancy comes from a failure of our approach 
in this regime, from other scattering mechanisms, or is 
an artifact of the deconvolution procedure which was not 
ruled out in Ref. In the present work we interpret 
the results using a spin 1/2 model which, as we believe, 
is valid at least when the renormalized coupling constant 
is small, C 1. This regime is usually referred to as 
"above the Kondo temperature" and its meaning in non- 
equilibrium situations will be explained in the discussion 
below. 

For vanishing magnetic fLeld the experiments by Poth- 
ier et alt! and Pierre et alB have shown that the distri- 
bution function depends only on the ratio e/elJ rather 
than on e/ksT. This observation is compatible with the 
fact that the Korringa inverse time, K, is proportional to 
eU at eU 3> fceT. Moreover, using the results in Sec. [II 
we would even get a distribution function which depends 
on the pair of dimensionless energies, e/eU, and, Eh /eU, 
only which qualitatively meets the experimental findings. 
This relation, however, relies on the assumption that the 
coupling constant is isotropic and energy independent. 



A. Inelastic relaxation rate 

Since the distribution functions both of electrons and 
spins are out of equilibrium, the inelastic relaxation can- 
not, in general, be discussed in terms of a single relax- 
ation rate. In the following we use two quantities to 
describe the inelastic relaxation. 

First, we consider the collision integral (0) in relaxation 
time approximation. We find the corresponding rate to 
be 



G 



— = / de>W(e'){l-f(e-e') + f(e + e>)} 

Trt J 



de'f(e-e')[W(e')-W(-e')], (19) 



with the spin-flip time defined by the exact relation fol- 
lowing from the sum rule ( p^ ) for the correlation function, 

C(u>), 



deW(t) = ^-^S{S + l){ P jf = —. (20) 



2h p 



This relaxation time equals the imaginary part of the 
electron self energy considering only interaction with lo- 
calized impurity spins. To discuss this expression at fi- 
nite magnetic fields we approximate W(e) by its weak 
coupling form to get 

— = ^{l-\[f{e-E H )-f{e + E H )] [V+-V-) 

(21) 

where (V+ — V-Yis the mean polarization of the localized 
spin, see Eqs. (14|) and (|l5|). For monotonous distribu- 
tion functions the additional term lead to a non-positive 
contribution and we have always l/r rt < 1 /r B f. At small 
magnetic fields the combination of occupation numbers 
and distribution functions in (|2l]) decreases and becomes 
of order (Eh /ell) 2 . In equilibrium we find this contribu- 
tion to be ~ tanh 2 (/3i?ff /2) that leads to an exact can- 
cellation of the spin-flip term, l/r rt = l/3r s f, for large 
magnetic fields, Eh 3> fc^T. This is in accordance to 
the fact that for large magnetic fields the spin-flip con- 
tribution is completely frozen out. Out of equilibrium, 
the contribution stemming from the distribution func- 
tions behaves similarly and cannot explain the increase 
of energy relaxation for small magnetic fields. 

Whereas in the collision operator (||) small energy 
transfer is canceled by the distribution functions the re- 
laxation rate in Eq. ( |1S| ) includes terms stemming from 
elastic scattering for e' = 0. To consider the pure inelas- 
tic part we define an inelastic scattering rate by 



1 



1 

Trt 



1 

T e \ 



1 

Trt 



K' 



de'W(e') 



(22) 



K' 



with some given cutoff, K' , not further specified. 

To understand the behavior of the inelastic relaxation 
time at small magnetic field one has to recall that the 
scattering processes at B = are almost elastic, a typical 
energy transfer being smaller or of the order of the Kor- 
ringa rate, K. Thus, subtraction of the elastic processes 
leads to a drastic decrease of the inelastic relaxation rate, 
or an increase of the inelastic relaxation time. For finite, 
small magnetic fields, however, the inelastic relaxation is 
of order K + Eh and therefore will lead to an increase 
of order (Eh/K) 2 . This term has to be compared with 



the monotonous decrease of order (Eh /elf) 2 and domi- 
nates in the regime p J <C 1 explaining the increase of the 
inelastic relaxation rate for small magnetic fields. 

A second way to discuss inelastic relaxation is to in- 
troduce the so-called energy relaxation time defined as 
an average energy transfer rate. It is defined as 



- = --_( dep(e)el{f}. 

Te e J 



(23) 



Here p(e) is the electron density of states, while I is a con- 
stant reference energy to normalize the particle energy. 
A natural scale for e at eU » fc^T is eU. Since eU <C D, 
where D is the electronic bandwidth, the electron den- 
sity of states can be regarded as energy-independent, and 
after some algebra we get 

~ = -~ JdeJduW(u))uf(e)[l-f(e-u)]. (24) 

One notices that this expression does not contain elastic 
scattering. 




FIG. 4. Inverse energy relaxation time l/r e (solid line) in 
arbitrary units for U = 0.3mV depending on magnetic field 
B measured in T. The dashed line shows the same quantity 
rescaled by 2/30 where all transferred energies ui in Eq. ( pl| ) 
are added positively. 



In Fig. |j we show the magnetic field dependence of 
the inverse energy relaxation time 1 /r e in arbitrary units 
for U = 0.3mV (solid line). We observe a monotonous 
decrease of the total energy transferred with increas- 
ing magnetic field. Considering Eq. ( |24| ) as the aver- 
age (p/e) J dedui ... of the energy transfer, one finds in 
the integrand, depending on the energy, both positive 
and negative contributions. Although the average j24| ) 
monotonously decreases when B increases, the positive 
and negative contributions in fact increase for small mag- 
netic fields Eh <C eU. For comparison we show by the 
dashed line the average ( |2^ ) where all transferred energies 
lo are added positively. This quantity, in contrast, show 
the non-monotonous behavior in dependence of magnetic 
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field. The physical relevance of this quantity, however, is 
not clear. This observation shows that the notion of a sin- 
gle energy relaxation time to describe energy relaxation is 
quite misleading. What happens is that non-equilibrium 
spins just redistribute the electron energy between dif- 
ferent energy regions, the extra energy being transferred 
due to diffusion of "hot" electrons in real space. 

In the experimental papers Refs. [l as m 

the theoretical ones based upon the i-matrix approach, 
Refs. H |10[ a contribution of two-particle scattering 
processes was also discussed. The results are expressed in 
terms of an effective interaction strength 7 which in our 
notation reads 7 = (pJ) 2 /r s f. In the experimental papers 
the results for B = and for finite magnetic fields were 
fitted assuming the same shape of the electron-electron 
interaction kernel oc 1/uj 2 . According to our findings, the 
behavior of the interaction kernel differs from ui~ 2 . Con- 
sequently, it is difficult to compare their and our results. 



B. Parameters obtained from experiment 

Now let us briefly discuss the parameters obtained from 
the fit of the experimental data. Similarly to Ref. where 
the first experimentstil with Cu were fitted, we determine 
the Kondo temperature and the impurity density (see 
Sec. IV) from the distribution function at B = 0. These 
values yield Tk ~ 0.4 K and c- lmp = 8 ppm, respectively. 
The procedure proposed in Ref. [12], based on the analy- 
sis of the saturated phase breaking time, would suggest 
a Kondo temperature below T « 0.1 K, see also Ref. |l3| . 
However, assuming such a small Tk in our formalism, 
we would disagree with the experimentally observed dis- 
tribution function at B = in Ref. |ll|- Here, we want 
to mention that we used the Kondo temperature in its 
simplest form based on the leading logarithmic approxi- 
mation. Corrections may still be large for this quantity 
and this disagreement could possibly be resolved using a 
refined version. 

The procedure based on the estimate of the deco- 
herence time, T<f,, leads to some other inconsistencies. 
Namely, an estimation of the phase coherence time with 
our parameters would lead to w T s f/2 « 0.06 ns while, 

In return, a fit of the de- 



according to Ref. 12 



1 ns. 



coherence time in Ref. 13 using the same r s f has lead to an 
impurity density of Ci mp ~ 0.15 ppm much lower than our 
estimated impurity density. This means that using the 
standard theory of weak localizationt3 the impurity spins 
should produce much higher decoherence rates than ob- 
served experimentally. The main problem in comparing 
these two quantities may be the following - the decoher- 
ence time is measured in equilibrium and it saturates in 
the case of Cu below the Kondo temperature if we assume 
that both and the energy relaxation are determined by 
the Kondo impurities. In that region, however, our treat- 
ment, as well as the standard theory of weak localization 
using the spins as random classical objects, become inac- 



curate. It is the aim of future work to develop a theory 
valid below Tk- 

In our treatment we used some assumptions which were 
not addressed so far. First, we used the same distribution 
function for the electrons with different spin projections 
and the standard diffusive form of the Boltzmann equa- 
tion. The latter assumption can be justified because the 
elastic relaxation timeoa Tj mp w 0.01 ps is much smaller 
than that for inelastic scattering. The usage of a single 
distribution function is certainly fulfilled when the spin- 
orbit relaxation rate, t" 1 , is larger than the rate t~7 of 
processes which tend to violate the electron spin symme- 
try. This assumption is supported by the results of the 
experimenttd yielding r so w 39 ps whereas r s f is of the 
order of few nanoseconds. 

We have disregarded some scattering mechanisms to 
focus on the magnetic dependence and to discuss the in- 
fluence of this sort of impurities only. Direct electron- 
electron interaction leads to a pronounced smearing of 
the distribution function at low energies. That can be 
important at rather large magnetic fields when the en- 
ergy relaxation mediated by magnetic impurities van- 
ishes. Since this additional effect at B = is much 
weaker than the contribution by magnetic impurities it 
cannot explain the deviations from experiment in the re- 
gion of small energy transfer. 

In the Appendix we determine the renormalized cou- 
pling constant J which replaces the bare coupling con- 
stant Jo- The calculations are performed within a RPA- 
like approximation. All higher order corrections are sup- 
posed to be taken into account by the renormalization of 
the coupling constant J. This procedure, as known, leads 
to the so-called overcounting problem which has been, 
considered both for the electronEZI and pseudo Fermiorcij 
self energy. In these papers an attempt to cure this prob- 
lem by setting the energy variable of the renormalized 
quantity to some special value has been made. No gen- 
eral method that would allow to avoid the overcounting 
has been suggested so far. The only message is that as 
long as the renormalization is weak and the renormal- 
ized vertices are changing slowly with energy one can 
use the "double dressed" vertices instead of the "single 
dressed" ones, using the language of Ref. |27|. Here, we 
simply adopt this "rule", however, having in mind that if 
the renormalized quantities get too much structure, i.e. in 
the Kondo regime, this approach may become inaccurate. 
Our renormalized vertices are similar to those obtained 
using the Hamanp-approximationcij in the i-matrix ap- 
proach by KeiterE2l,_which is a generalization of Suhl's 



dispersion approached to finite magnetic fields. We fur- 
ther neglected the non spin-flip contribution which may 
be put into the free Hamiltonian and does not influence 
energy relaxation. 

In the calculation of the correlation functions (||) 
and ([l0]) we neglected the Knight shift as well as the 
frequency and Zeeman energy renormalizations. The 
Knight shift just leads to a simple addition to the "not 
really known" Zeeman splitting of the impurity spin and 
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can be disregarded. The other terms can be deter- 
mined by the Kramers-Kronig relation from the Korringa 
widths (JlJ) and ([l^). They vanish for e, Ejj = and give 
just a renormalization of higher order in (pJ) 2 which is 
usually disregarded well above the Kondo temperature. 
However, below the Kondo temperature these corrections 
would lead to an anomalous behavior and one has to go 
beyond the lowest order in the memory function even 
if one uses already renormalized coupling constants, see 
Ref. for a discussion. 



C. Kondo Temperature 

The regime where all our calculations are actually valid 
is defined here as the regime where the lowest ord er in 
t he logarithmic corrections, such as in Eqs. (54) and 
(B5), are still small compared to 1. This is what we ac- 
tually mean when speaking about a "non-equilibrium" 
regime or about a situation "above the Kondo tempera- 
ture" . The Kondo temperature is usually defined by the 
electron temperature where perturbation theory starts to 
fail. The terms of the perturbative expansion explicitly 
depend on the electron distribution function. As a result 
even when the bath temperature is reduced below Tk the 
perturbative approach can still be applied provided that 
the deviation from equilibrium is strong enough. 

In Refs. H ||, and [h] the use of an effective spatially de- 
pendent Kondo temperature T£(x) = T][ x /(ef/)^"^ 
was proposed for electrons energetically near the Fermi 
point e = — ell/ 2 using the free solution (|l|) as a basis of 
the renormalization. For electrons near e = eU/2 one has 
simply to replace x — > 1 — x. In these terms the effective 
Kondo temperature in the non-equilibrium situation is 
always smaller than the equilibrium Kondo temperature. 

For an arbitrary distribution function it is not as 
straightforward to calculate the renormalized quantities 
as it is for the free solution (|l|). Equivalently to the 
notion of a renormalized Kondo temperature we pro- 
pose to use a local temperature T*(x) describing the 
non-equilibrium situation. The Kondo temperature then 
equals its bulk equilibrium value and the local temper- 
ature is a spatially dependent functional of the distri- 
bution function. For strong electron-electron interac- 
tion or equivalently in the middle of a very long wire 
T*(x) equals-the|-apalytically determined "hot electron" 
temperaturet3'tZHl3 independent of the scattering mech- 
anism. Further, this temperature fulfills the boundary 
conditions T*(x = 0,1)= Tbath- Therefore, at the ends 
of the wire and for temperatures explored experimentally 
we are below Kondo temperature and our approach is 
inadequate. Here, we assume that the boundary regions 
where our approach fails are narrow and do not influ- 
ence the distribution function in the middle of the wire. 
We believe that this assumption is fulfilled rfor long wires 
that were studied in the experiments.EJOlI.3 Further, it 
is obvious that our approach becomes more accurate for 



larger applied voltages when the distribution function at 
B = is a function of e/ell only. This is the reason why 
we focus more on U = 0.3 mV. It turns out (see below 
for details) that for U = 0.1 mV and smaller voltages our 
theoretical approach reaches its limit of validity within 
the experimentally explored temperature range. 

To be more specific we give an analytic expression for 
the effective temperature in the case of weak smearing. 
The distribution function then equals the free solution 
(|]) and the smearing is just included in an effective tem- 
perature T c ff characterizing the energetic width of the 
smeared steps. (Here, the Boltzmann constant ks is set 
to one.) For simplicity we restrict ourselves to B = but 
the generalization is straightf orw ard. A s a basis of the 
renormalization we use Eqs. (BS) and (|i?5| ). Following 
the usual poor man's scaling procedure we get 



pj{e) 



pJo 



1 - pJoff(e) 



(25) 



whereby we adiabatically integrated out all high energy 
contributions from —D to e — D and D to e + D letting 
D tend to zero. This procedure at finite temperature 
or finite effective temperature leads to an effective lower 
energy cutoff which is either e, eU/2, or T c ff and follows 
immediately from 



fl (e)w(l-i)ln 
+ x In 



D 



\e-eU/2\ + T eS 
D 



e + eU/2\+T ctf 



(26) 



Inserted into ( pq) we get for the renormalized coupling 
constant 



pj(e) = ln 



with 



{je - eU/2\ + Teff) 1 -^^ + eU/2\ + T oS ) x 
T K 



T, 



K 



De 



-1/pJo 



(27) 



(28) 



For sufficiently weak smearing T c ff <C eU mostly elec- 
trons near the Fermi points e ±eU/2 contribute. We 
may define a local temperature for electrons depending 
on the Fermi point. For electrons at e ~ eU/ 2 we get 
T+(x) = T^ x (eU) x . The other local temperature for 
electrons at e « -eU/2 reads TZ(x) = T^el/) 1 '* . 
Here, of course, the effective smearing T c g depends on x 
and coincides at the ends of the wire (x = 0, 1) with the 
bath temperature Tbath- Since electrons of both Fermi 
points contribute to inelastic processes one has to average 
the effective temperatures. The simplest average which 
yields the correct expression in the middle and the ends of 
the wire reads T*(x) = (1 -x)T^ x {eU) x +xT* s {eU) 1 ^, 
because at x = 0, 1 only electrons at one Fermi point con- 
tribute. 







In the middle of the wire we get simply the require- 
ment T c geU 3> T\ for the calculations to be valid. If 
the smearing becomes stronger and T a is of the order of 
eU or even larger we get from Eq. (|27| ) the requirement 
T* = T cff » T K . 

On the other hand, we may define an effective Kondo 
temperature by the failure of the perturbation theory. 
For electrons near the lower Fermi point e w —eU/2 we 
write 



pj{e) 




\eU\( 1 ~ x ' ) / x 



T, 



l/x 



K 



(29) 



and find the effective Kondo temperature to be T^(x) — 
t k X /{eU)^- x ^ x as proposed in Ref. 0. The Kondo 
temperature for electrons near the upper Fermi point fol- 
lows analogously. To our opinion the introduction of a 
local temperature is much more intuitive than a modi- 
fication of the Kondo temperature. It seems to us that 
the latter concept can sometimes be misleading. 



To show that magnetic impurities are indeed present in 
copper wires one has to describe all the available experi- 
ments using the same set of parameters. Since there are 
strong deviations to the magnetoresistance experiments 
this is still an open question. 
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APPENDIX A: DERIVATION OF THE 
COLLISION INTEGRAL 



D. Conclusion 

In this article we studied energy relaxation mediated 
by magnetic impurities above Kondo temperature. We 
have shown that the effective electron-electron interac- 
tion is already included in the dissipative nature of the 
spin system. The finite line width in non-equilibrium 
proportional to the applied voltage lead to strong en- 
ergy relaxation even for vanishing magnetic field where 
the spin states are degenerate. We succeeded to derive a 
non-perturbative description valid for arbitrary magnetic 
fields where the perturbative t-matrix approach fails due 
to divergences in the two electron scattering rate. 

We characterized the efficiency of the energy relax- 
ation by the slope of the distribution function at e « 0. 
The magnetic field dependence of this quantity shown 
in Fig. shows a non-monotonous behavior. For small 
magnetic field the energy relaxation strength is increased 
and for larger magnetic fields it decreases again. These 
features meet qualitatively recent experimental data on 
energy relaxation in mesoscopic copper wires. Although 
the detailed comparison of the energy distribution func- 
tions lacks to fit quantitatively all energy regions, this 
gives a strong indication that magnetic impurities are 
indeed responsible for energy relaxation in copper wires. 

Several things remain to be done in order to achieve a 
quantitative description for the energy relaxation exper- 
iment in the broad regions of temperature and magnetic 
field. To go below the Kondo temperature the overcount- 
ing problem has to be addressed. High-order corrections 
in impurity density might require a more complicated 
kinetic equation including quantum effects. Finally, for 
larger magnetic fields and/or low concentration of the 
localized spins other relaxation mechanisms have to be 
accounted for. 



We start with the 'bare' Hamiltonian H = Hq + Hi 
where the free Hamiltonian is given in Eq. (Q). In terms 
of pseudo Fermions it reads: 

h = J2 ^°L C ^ + E E p4 a P • ( A1 ) 

kcr /3 

Here the creation/annihilation operators Ct ,C k de- 
scribe electrons while the operators ttgjOg describe 
pseudo Fermions. The pseudo Fermion states are speci- 
fied by (3 = ± having the energies E± — ^fgfiBB/2. The 
interaction Hamiltonian 

Hi = J Q S ' S v'° C l>a> C kcr 

kk' aa' 

is the s — d exchange Hamiltonian with the bare coupling 
Jo- Expressing spin operators through pseudo Fermions 
one can rewrite it as 

kk' 

(summation over repeated indices is expected). In addi- 
tion, we have to take into account the operator constraint 
a} + a + + a}_a_ = 1 which can be formally done by a pxo- 
jection with the help of a complex chemical potential^. 

Linear in the impurity density Ci mp the angular aver- 
aged collision integral fcp,the classical Boltzmann equa- 
tion can be expressed astSa 

/{/} = ^{/(c)S > (e) + [1 - /(e)]S<(e)| (A3) 

with S >/< (e) = J^a S >/< (kcr,e)/2 where e = e k(T is the 
spin averaged self energy assumed to be independent of 
the angular momentum, /(e) is the angular averaged 
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distribution function for electrons of energy e where we 
suppressed the spatial dependence for convenience. Since 
the self energy is proportional to the impurity density we 
already replaced the electron Green's functions by their 
unperturbed form and integrated over frequency to get 
the classical form of the Boltzmann equation. 



FIG. 5. Second order self energy graph for the electron 
Green's function 



The electron self energy in lowest nonvanishing order, 
depicted in Fig. ^[ is given by a pseudo Fermion bubble 
and an electron or hole line in between 

S>(ka,e) = c imp j5^5E / ?^>(kV, f - W ') 



electron-hole bubble are depicted in Fig. [?]. Considering 
only these two topologically different graphs in the self 
energy we obtain the i-matrix expression for the collision 
integral derived in Refs. however, with bare coupling 
constants. The electron-hole pairs do not only renor- 
malize the pseudo Fermion propagator, which has been 
considered in Ref. |2^, but also include bubbles connect- 
ing the upper and lower pseudo Fermion line. The latter 
graphs are so palled crossed diagrams and lead to ver- 
tex corrections! 2 !] not included in a simple non-crossing 
approach. 

All other corrections lead to higher orders in Jo for a 
single electron-hole bubble or outer electron line. Above 
the Kondo temperature it is conventionally accepted that 
one may write these corrections as a renormalization of 
the corresponding vertices and that in the leading log- 
arithmic—order the vertices are renormalized indepen- 
dentlyBL3. As a result, two dressed vertices instead of 
one are used. This procedure has been explicitly proven 
in the leading logarithmic approximation both for the 
electrorO and pseudo Fermiono self energy. An alterna- 
tive derivation of the independent vertex renormalization 
has been given in Ref. [?]. 



(A4) 



where summation over internal spins is implied. The self 
energy S, is given by the obvious change of the < and > 
signs in (A4). Since all corrections in the electron Green's 
functions are of higher orders in c; mp they can always be 
used in the free form G^kcr, u>) = 2-Kif(uj)5(uj — £k<r) 
and GqQ/uj,uj) = —2ni{l — f(oj)}S(ui — e k(T ), respec- 
tively. The pseudo Fermion Green's function in low- 
est order reads Q^-(f3,u>) — 2niT :> pS(Lu — Ep) and has to 
be renormalized in the follqsdng. Note that when us- 
ing the Abrikosov technique^, the occupation numbers 
V± acquire an additional chemical potential exp(— iX) 
and therefore Qq(/3,uj) = —2-Ki5(oj — Ep), see Ref. |3^. 
In contrast to Ref. [32] we assume that the occupation 
probabilities V± are determined by the non-equilibrium 
electron distribution. Therefore they have to be found 
selfconsistently. 



FIG. 6. Fourth order graphs which are of second order in 
impurity density 



For further proceeding we have to classify the appear- 
ing graphs. Graphs including two pseudo fermion bub- 
bles like shown in Fig. |^ are of second order in the im- 
purity density and are therefore neglected. It is obvious 
that we may include an arbitrary number of additional 
electron-hole bubbles where each one is at least of or- 
der Jg . The two lowest order graphs with one additional 



FIG. 7. Fourth order self energy terms including one addi- 
tional electron-hole bubble. 



Assuming this independent renormalization of vertices 
to be correct also for more involved graphs and that this 
renormalization depends only on the electron energies, 
one can equivalently start with a renormalized Hamilto- 
nian restricting the appearing graphs to simple electron- 
hole bubbles of the order of J 2 in the renormalized matrix 



elements Jf,f . 



Neglecting the energy dependence of one obtains 



E > (kcr,w) : 



„. ffi'fi J/3/3' V- 



duj' du >f , , , 



<e > (/3 / ,^ + ^)eo < (A^)) ch 



(A5) 



where the average (. . .) c h means the dressing of the 
pseudo Fermion bubble with all possible electron-hole 
pairs using renormalized coupling constants. Rephras- 
ing the pseudo Fermions in terms of spin operators and 
assuming isotropic coupling J we may write 



4mp , 

~4p~ 



{pjf / duj'C{u')[\ - f(w - uj')] (A6) 



with the time-dependent correlation function C (t) given 
by Eqs. (||) and (^). The self energy E < is given by the 
replacement / — » 1 — / and C(ui) — > C(— uj) which is ob- 
vious since C (t) is an autocorrelation function. Inserting 
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this expression in (A3) we find the collision integral in the 
desired form . The spin-spin correlation functions a*e 
calculated then using a projection operator techniqucEil. 
The actual procedure follows the lines outlined in Rcf. |36| 
for the case of two-level systems. 



APPENDIX B: VERTEX RENORMALIZATION 

We proceed to determine the renormalized coupling 
constants. They can be derived—ip different ways: 
by a pool, man's scaling procedureO, Suhl's t-matrix 
approach!!!!, following. Abrikosov and solving a general- 
ized vertex equatior£3, or considering lower order correc- 
tions and assuming similar resummation procedure as in 
equilibrium. According to Refs. ^7] and ^ the renormal- 
ized vertex reads 



r(eV, oj'(3'\ea, w&) = T Q + T e (e + u) + T h (e - 



where 



To = J 



■/3'P 



(Bl) 



(B2) 



is the bare vertex of the s — d exchange Hamiltonian. a, (3 
are the incoming electron and pseudo Fermion spins and 
e, uj are the incoming electron and pseudo Fermion ener- 
gies. The primed quantities are the outgoing spins and 
energies, respectively. Energy conservation is fulfilled at 
each vertex meaning e' = e + uj — uj' . The electron and 
hole vertex parts can be assumed to depend on a sin- 
gle energy variable.ES The electron vertex r e depends on 
the sum of incoming electron and pseudo Fermion ener- 
gies e + oj, and the hole vertex part T/,. on the difference 
of the incoming electron and outgoing pseudo Fermion 
energies e — uj' . In the lowest order in Jo the retarded 
quantities read 



2tt 



G r (qK,e T 0^(7,0 



Here the upper sign stands for the electron vertex and the 
lower sign for the hole vertex, respectively. The auxiliary 
function 



9(e) 



_, f(e')-l/2 _ 
- D e — e' + iS 



dc 



(B3) 



leads in equilibrium to the usual logarithmic corrections. 

Considering the lowest order vertex corrections to the 
self energy in Fig. || and assuming that the pseudo 
Fermion energies are fixed to resonance, we get 



Jo 



T ± = y i [ l + ^[g(e±E H )+g(e)} 
for the spin-flip component proportional to and 



(B4) 



T z ± = (J /4) [1 + pJ Q g(e ± E H )\ 



(B5) 



for the component proportional to S z applying to a spin 
up/down electron. This means that the leading renor- 
malization of the 'up' electron is given by the 'down' elec- 
trons and vice versa. This anisotropic decomposition of 
the renormalized coupling constants has also been found 
by Keitera. Setting the pseudo Fermion energy to res- 
onance is correct only for the lowest order correction in 
the electron self energy in Fig. However, considering 
corrections to the graphs in Fig. ^ we find that there 
the pseudo Fermion energies are weighted with a diver- 
gent term oc (uj — Ep)~ 2 justifying the use of the reso- 
nance energy uj w Ep in the renormalization. Assuming 
that this is also justified for more involved graphs we use 
these corrections as a basis of our renormalization. To 
write the renormalized Hamiltonian in terms of effective 
anisotropic and energy dependent coupling constants we 
reintroduce the spin operators and define J±(e) = 4Fj_(e) 
and J ± (e) = 2r ± (e). 



FIG. 8. Lowest order vertex renormalization 



To define a 'Kondo scale' we consider the renormal- 
ized coupling constants as functions of the incoming-elec- 
tron energy and use the Hamann approximationE3 for 
the renormalized quantities. This approximation yields 
the same leading-logarithm expansion as the p oor m an's 
scal ing procedure based on the expressions (B4) and 
(B5), however, taking care of the unitarity condition in 
each order. In this approximation the S z coupling con- 
stant reads 



J|(e)/J = { |1 - (ir P J ) 2 S(S+ l)/4- P J g(eTE H )\ 2 
+ (7rpJ ) 2 S(S+l)y 1/2 (B6) 



in accordance with a high temperature expansion of 
KeitercJ. Analogously the spin-flip process renormaliza- 
tion reads 



J 



z (e) /Jo = {|1 - (ir P J ) 2 S(S+ l)/4 - pJ \g(e) 
+g(e ± E H )}/2\ 2 + (tt P J ) 2 S(S + 1)}~ V ' 



(B7) 



These quantities coincide for B = with those used in 
Refs. @|| and cover the results obtained by the poor 
man's scaling procedure, cf. discussion about Kondo tem- 
perature in Sec. 0. Note that the renormalized coupling 
constants remain finite also below the Kondo tempera- 
ture where their meaning is questionable. 
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APPENDIX C: CALCULATION OF THE 
CORRELATION FUNCTION 

We briefly present the method used to derive the re- 
sults in Sec. [II. Our main interest is the spin-spin 
correlation functions (^). Since C(t) is an autocorre- 
lator its Fourier transform is given by the expression 
C(u) = C(-iu) + 5) + C{-iu) - 5) where C{z) the 
Laplace transform of the correlator. Using projection 
operators P Z X = S z (X S z ) / (S z S z ) for C z and P ± X = 
S (X S T ) I (S^ S T ) for the C± component, one can de- 
rive a formally exact integro-differential equationtiJ 

Ca{t) = ^aC a {t) - f d U( f> a (t ~ u)C a {u) (CI) 

Jo 

with the solution in terms of the Laplace transform 

C o (i = 0) 



C a (z) 



z - <f> a + 4> a (z) 



(C2) 



where a = z,±. Here, $ z = (S* S z ) / (S z S z ) = and 
$± = (S ± S^)/(S ± S^) = TiE H leads to the free propa- 
gation of the correlators, where Eh includes the Knight 
shift neglected throughout the article. The averages are 
calculated with the proper steady state density defined by 
the stationary solution of the Boltzmann equation. The 
memory kernel <p a it) for the correlation function plays a 
similar role as the self energy for the Green's function. 
For the C± term we find 



(C3) 



Here, the index r in S^(t) indicates that the dynamics of 
the spin operator is reduced by the projection. It is de- 
termined by the expression S^(t) = exp[iL(l — P ± )t]5 ± 
with the Liouville operator L acting as LX = [H,X]/U. 
The memory ker nel for the C z correlation function is 
simply given by (C3) with the replacement ±, =p — > z. 
Though formally exact, the above equations do not al- 
low a simple calculation of the correlators. Here, we ex- 
pand the kernel up to second order in the renormalizcd 
coupling J. Since the dynamics of the expanded kernel 
function is oscillatory the Fourier transformed correlation 
function has always the simple form 



C a (L0) = 



2C a (t = 0) Re^ a (u;) 



[u - i$ a + Ini0 a (w)] 2 + [Re0 Q (w)] s 



(C4) 



with a = z, ±. Similar to Green's functions, this re- 
flects the fact that in the steady state the retarded and 
advanced self energies are complex conjugate, leading 
to a similar structure of the ">" or "<" Green's func- 
tions. Note that the roles of the imaginary and real 
parts of the memory kernel are opposite to those of the 
Green's functions because of different set of definitions. 
Further, we define Ke(f> a (uj) = Re{(/> a (— iui + 8)} and 



the imaginary part Iia<j) a (ui) follows from a Kramers- 
Kronig relation. With C z (t = 0) = (S Z S Z ) = 1/4 and 
C±(t = 0) = (S ± S T ) = V± and further neglecting the 
imaginary parts in the denominators which lead to a fre- 
quency and Zeeman energy renormalization, we find the 
expressions (0) and (fill) for the correlation functions. 



APPENDIX D: DETAILED THEORETICAL 
DESCRIPTION 

Starting with anisotropic and energy dependent cou- 
pling constants derived in Appendix [b| the interaction 
Hamiltonian in terms of impurity spin operators reads 



kk' 



-S' 



[j^(e feT )4 T C7 fcT -Ji(e fci )4 i ^ 



(Dl) 



Considering the electron self energy in the Boltzmann 
equation the scattering rate W in the collision integral (||) 
depends on both, incoming electron energy e and trans- 
ferred energy u> 

W(e,u) = (c impP /8h) {j-{e)J+(e')C+{uo) 

+J + (e)J-(e')C-(uj) (D2) 
+ [J+(e)J+(e') + Jt(e)Jt(e')] C z (uj)} 

with the outgoing electron energy e' = e — to, cf. Ap- 
pendix |a| and Eqs. (||) and (||). The calculation of the 
correlation functio ns f ollows the lines in Appendix K2 and 
the general form (C4) leading the Eqs. (|j) and (|lC) re- 
main the same. The changes are only in the damping 

v z (io) = irp 2 [T + C+(uJ - E H ) + V-C-{u + E H )\ (D3) 



i/±(w) = (tt/4) P 2 [C z (w T E h ) + t T (u)/V± 



(D4) 



with the auxiliary functions 

U")=jde' [J z + {e')j z + {uj + e') + J*_(e')J z _(u, + e')] 

xf(e')[l-f(uj + e')} (D5) 
( ± (Lj) = [de'J*(e')J ± (u + e')f(e')[l - /(w + e')] . (D6) 



The corresponding master equation ( |i4[ ) for the occu- 
pation probabilities V± remains the same but the rate 
changes to 



r± = 



p 



4frP± 



duj Q± {—uj)C± (to) . 



(D7) 



We used these formulae for the numerical determination 
of the distribution function in Sec. |iv| , however, the main 
features do not alter even if one uses the simplified de- 
scription in Sec. III. 
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